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Inelastic scattering induces dephasing in mesoscopic systems. An analysis of previous models to 
simulate inelastic scattering in such systems is presented and also a relatively new model based on 
wave attenuation is introduced. The problem of Aharonov-Bohm(AB) oscillations in conductance 
of a mesoscopic ring is studied. We have shown that conductance is symmetric under flux reversal 
and visibility of AB oscillations decay to zero as function of the incoherence parameter, signalling 
dephasing. Further wave attenuation is applied to a fundamental problem in quantum mechanics, 
i.e., the conditional(reflection/transmission) times spent in a given region of space by a quantum 
particle before scattering off from that region. 
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We present a study of two different problems in mesoscopic physics which are of current interest namely, dephasing 
and the other conditional times. This work considers these two phenomena not from a microscopic point of view 
but through a phenomenological model which captures the essence of these quite well. This model is known as wave 
attenuation. It essentially involves damping the wave function in the region of interest, so as to derive essential physics. 
In the context of dephasing, we consider an Aharonov-Bohm interferometer, and show that the wave attenuation model 
is better than its counterpart the optical potential model used in this context. Also, in case of double heterostructures 
we employ this technique to clock the time a particle takes to traverse a local region of interest in the given system. 
Here also the wave attenuation model is more simpler to deal with than other models. 



The process of dephasing or decoherence leads to the diminishing of quantum effects or loss of quantum mechanical 
interference effects. Dephasing occurs due to interaction of an electron (interfering entity) with it's environmental 
degrees of freedom (which are not measured in the interference experiment) P, • The AB oscillations are one of 
the prime examples for analyzing how quantum interference effects are affected by dephasing. These oscillations are 
similar to the fringes seen in an Young's double slit experiment apart from the presence of the magnetic flux. In 
the Young's double slit interferometer the intensity is given by I = l^ 2 = \ipi\ 2 + \ip2\ 2 + 2Sft( , 0* V , 2 e *^)j the part 
25R(-05'V'2 eI ^) represents the interference term. Here ip\ and ip2 ar e the complex wave amplitudes across the upper 
and lower arms of the interferometer and <f> is phase difference between these two wave amplitudes. If there would 
have been no phase relationship between the waves then the average intensity will be < I >= \ipi\ 2 + \ip2\ 2 ■ Complete 
dephasing is indicated by extinction of these interference terms. Thus, Dephasing can be defined as a phenomena by 
which quantum mechanical system behave as though they are described by classical probability theory. 



In absence of a complete microscopic theory as to how inelastic scattering affects dephasing, models are useful. 
There are different ways to model dephasing in mesoscopic systems. An interesting method is to attach a voltage 
probeQ to the sample as in inset of Figure 1 (Buttiker's model). In this model, an electron captured by a voltage 
probe is re-injected back with an uncorrelated phase leading to irreversible loss of phase memory. This model has 
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built in current conservation and Onsager's symmetry relations are obeyed but it does suffer from a major demerit in 
that it only describes localized dephasing. 

The optical potential model provides another method of introducing dephasing in these mesoscopic samples. The 
optical potential was first introduced to explain the inelastic crossection in case neutron scattering. In this model 
a spatially uniform optical potential(— iVA is added to the Hamiltonian£|, making it non-hermitian. This leads 
to removal of particles from phase coherent motion, while absorption is real in the case of photons there is no 
absorption for electrons. The absorption of electrons is reinterpreted as scattering into other energy levels and 
therefore proper re-injection is necessary. Zohta and Ezawa 4] interpreted that the total transmission is defined after 
re-injection as the sum of two contributions one due to the coherent part and the other due to the incoherent part, 
i.e., T tot — + Tj nco / l . The incoherent part is calculated as T inco h — T T ^ T A, herein T r and T\ are the probabilities 
for right and left transmission from the region of inelastic scattering and A is the absorbed part which is given by 
A = 1 — T co h — R co h- This model unlike Buttiker's model describes dephasing that occurs throughout the system but 
still it has some major drawbacks in that it violates Onsager's two terminal symmetry relations 6J and in the limit of 
strong absorption leads to perfect reflection. 





FIG. 1: Aharonov-Bohm ring geometry. Inset shows the three probe model. 



Thus, there is need of a model which is free from the shortcomings of the above two models. Brouwer and 
BeenakkerQ have developed a simple method, by mapping the three probe Buttiker's method into a two terminal 
geometry, this is done by eliminating the transmission coefficients which explicitly depend on the third probe by 
means of unitarity of the S-matrix. They consider a three terminal geometry in which one of the probes is used as a 
voltage probe in absence of magnetic flux (see inset of Figure 1). A current I = I± = —I 2 flows from source to drain. 
In this model, a fictitious third lead connects the ring to a reservoir at chemical potential ^ in such a way that no 
current is drawn (1$ = 0). The 3X3 S-matrix of the entire system can be written as- 



S 



r\x tia *i3 

t21 ?~22 ^23 
t31 £32 ?~33 



Application of the relat ions 0,0)13 - 
two probe conductance G = — 3- 



I p = J2 q G P q[[ip—fJ.q},p — 1, 2, 3 and G pq = (2e 2 /h)T pq yields the (dimensionless) 



G = T 21 



T31 + T32 



where T pq — T p ^ q , is the transmission from q th to p th lead (\t p 



(1) 



,| 2 ), on elimination of the transmission coefficients 



in Eq. 1, which involve the voltage probe using the unitarity of the S - Matrix leads |3 to- 



G = T 21 
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Now all the above coefficients are built from the 2X2 S matrix- 

s' = ( : n ^ ) 

V hi r 22 J 

which represents the S Matrix of the absorbing system (non-hermitian) [i| . Thus re-injection has been reformulated 
as in Eq. 2. The first term in Eq. 2 represents the conductance contribution from the phase coherent part. The 
second term accounts for electrons that are re-injected from the phase breaking reservoir, thereby ensuring particle 
conservation in the voltage probe model. Also, Eq. 2 restores Onsager's two terminal symmetry relation for the optical 
potential model. 



B. Wave attenuation as a model for dephasing 



In-spite of the fact that a major problem associated with the optical potential has been cured there still remains the 
problem that in the strong absorption limit it leads to perfect reflection and absorption without reflection(spurious 
scattering) is not possible |lClllllll2| . To overcome this problem, we instead of making the Hamiltonian non-hermitian, 
add an exponential factor e~ al to the S-Matrix of the system. This is the model of wave attenuation. This describes 
dephasing which occurs throughout the system and removes the short comings of the optical potential model. 

The wave attenuation model is not new it has earlier been dealt with in the context of f-D localization^! ■ I 11 
the AB ring geometry considered here, wave attenuation is inserted by the factor e~ alt (or e~"' 2 ) in the complex 
free propagator amplitudes, every time we traverse 8] the upper (or lower) arms of the ring (see Figure 1). We have 
calculated the relevant transmission and reflection coefficients by using the S- matrix method along with the quantum 
wave guide theory for a single channel case. In this model, average absorption per unit length is given by 2a. With 
this method we show that the calculated conductance (G) in Eq. 2 is symmetric under the flux reversal as required. 
The visibility of the AB oscillations rapidly decay as a function of a, indicating dephasing. Hence forth we refer to a 
as an incoherence parameter. Increasing a corresponds to increasing dephasing processes in the system or increase in 
temperature. 

In Figure I, the length of the upper arm is l\ and that of lower arm is 1%. The total circumference of the loop 
is L = li + 1%. The loop is connected to two current leads. The couplers (triangles) in Figure 1 which connect the 
leads and the loop are described by a scattering matrix S. The S matrix for the left coupler yields the amplitudes 
0\ = (a', emanating from the coupler in terms of the incident waves I\ = (a, (3%, 71), and for the right coupler 

yields the amplitudes O2 — (8', /%, 7o) emanating from the coupler in terms of the incident waves I2 = (8, (32, 72)- The 
S-matrix for either of the couplers |l3( is given by- 

( -{a + b) yfe Vi\ 
S=[ yfi a b 
V y/Z b a J 

with a — (1 — 2e) — 1) and b = — 2e) + 1). e plays the role of a coupling parameter. The maximum 

coupling between reservoir and loop is e = ^, and for e = 0, the coupler completely disconnects the loop from the 
reservoir. 

The waves incident into the branches of the loop are related by the S Matrices for upper branch by- 

( Pi\ = ( e lkh e- ah e^ \ { P[\ 

{fc) ~ y e ikh e -ai le ^x J\&) 

and for lower branch- 




These S matrices of course are not unitary S(a)S{a)^ 7^ 1 but they obey the relation S(a)S(— ay = 1. The same 
relation is also obeyed by the S Matrix of the system in presence of imaginary potential. Here kl\ and H2 are the 
phase increments of the wave function in absence of flux. ^ and ^ are the phase shifts due to flux in the upper and 
lower branches. Clearly, _|_ ®h — where $ is the flux piercing the loop and $0 is the flux quantum^. The 

transmission and reflection coefficients in Eq. 2 are given as follows- T%i = \^\ 2 , Ru — l^| 2 ; R22 = ItI 2 ; T12 = lifl 2 
wherein 5' , <5, a' , a are as depicted in Figure 1 . 
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After calculating the required reflection and transmission coefficients wc sec that the coherent transmission T21 is 
not symmetric under the flux reversal however proper re-injection of carriers by Eq. 2 for the total conductance G 
plotted shows that the Onsager's symmetry relations are restored, i.e., G is symmetric under flux reversal. On can 
notice from this figure that amplitude of AB Oscillations decrease with increase in incoherence parameter a. All 
parameters used in the following figures are in their dimensionless form. 
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FIG. 2: Reduction of Aharonov-Bohm Oscillations in presence of incoherence parameter a. kL = n, h/L = 0.75, I2/L = 0.25, 
and coupling e = 0.5. 



In Fig. 3 we plot visibility (V) as a function of incoherence parameter a. Visibility is of course defined as- 

Gjxiax G m i n 



V 



Gmax ~\~ (~* mi 



(3) 




FIG. 3: Visibility for the same physical parameters as in Figure 2, coupling parameter e = 0.5. In the inset the harmonics have 
been plotted for the same physical and coupling parameters. 

The plot shows that with increase in the value of the parameter a the visibility exponentially falls off, reaching a 
point where it becomes zero corresponding to the disappearance of quantum interference effects, i.e., total dephasing in 
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the system. In the inset of Figure 3 we have plotted the first few Fourierp3 harmonics o„ (wherein n = 1 to 4) of G($) 
as a function of a. The harmonics are calculated from the formula-a n = J Q Gcos(n<f))d<f). The harmonics exponentially 
fall off with increasing a with exponent increasing as we go from 1 to 4. The n th order harmonic corresponds to the 
contribution from electronic paths which encircle the flux n times. The harmonics can sometimes show non monotonic 
behavior depending on the physical parameters|l5j|. however, the visibility is a monotonic function of the incoherence 
parameter a. 

Thus we have seen that wave attenuation can be effective in modeling absorption induced dephasing in mesoscopic 
systems. In the next section we will see how this method is extended to measure the time a quantum particle takes 
to traverse a specified region of scattering. 



III. TIME IN QUANTUM MECHANICS 



One of the most important problem in quantum mechanics is to calculate the time spent by a particle in a given 
region of space before scattering off from that region. The problem is essentially due to the fact that there is no 
hermitian operator to calculate this time in quantum mechanics ^| 0, [^. The prospect of nanoscale electronic 
devices has in recent years brought new urgency to this problem as this is directly related to the maximum attainable 
speed of such devices. When it comes to tunneling time or time in general there is lot of ongoing controversy. In some 
formulations this time leads to a real quantity and in others to a complex quantity jl9j. In certain cases tunneling time 
is considered to be ill-defined or quantum mechanics does not allow us to discuss this time^l |2(], H]]. Furthermore 
sometimes it is maintained that tunneling through a barrier takes zero time|22|. Recently, Anantha Ramakrishna and 
Kumar (AK)|23j have proposed the non unitary Optical potential as a clock to calculate the sojourn times without 
the clock affecting it. In this paper we examine another non-unitary clock, i.e., wave attenuation to calculate the 
conditional lengths, i.e., the total effective distance traveled by a particle in the region of interest. This conditional 
length on appropriate division by the speed of the particle in the region of interest will give us the conditional time. 

The criteria, any result for the time spent by a quantum particle in a given region of scattering should satisfy 
are that-(l) It should be real, (2) It should add up for non overlapping regions, (3) It should be causally related 
to the interval of space, and (4) tend to the correct classical limits. It is shown explicitly by AK that all clock 
mechanisms involve spurious scattering or very clock mechanism affects the sojourn times to be clocked finitely even 
as the perturbation due to clock potential is infinitcsimally small. All of these clocks involve spurious scattering as the 
perturbation due to clock mechanism couples to the Hamiltonian. This raises the question "Can quantum mechanical 
sojourn time be clocked with clock affecting it?" . In this paper we introduce such a method in which perturbation is 
not introduced in the Hamiltonian but in the S-Matrix of the system. In this case scattering is treated analogously 
with the Fabry-Perot interferometer. The scheme is illustrated in Figure 4. 



A. Wave attenuation to measure conditional times 

In the wave attenuation method^. Il0|. we damp the wave function by adding an exponential factor (e - "') every 
time we traverse the length of interest, here 2a represents the attenuation per unit length. This method is better than 
the optical potential model as it does not suffer from spurious scattering's^ 0, The corrections introduced in 
case of Optical potential model to take care of spurious scattering's will become manifestly difficult when we calculate 
the traversal times for a superlattice involving numerous scatterer's. Thus our method of wave attenuation scores 
over the optical potential model. In the presence of wave attenuation a wave attenuates exponentially and thus the 
transmission (or reflection) coefficient becomes exponential with the length endured in presence of the attenuator and 
this acts as a natural counter for the sojourn lengths. Following the procedure of AK we calculate the traversal and 
reflection lengths and times in a given region of interest (in particular between two scatterer's as in Figure 4. 

The amplitude for transmission and reflection can be calculated by summing^ the different paths as in Figure 4. 
The scatterer S± in Figure 1 has as its elements r 1 ,r' 11 t 1 and t\. r 1 is the reflection amplitude when a particle is 
reflected from the left side of the barrier while r( is the reflection amplitude when a particle is reflected from the right 
side of the barrier. t 1 and t[ are the amplitudes for transmission when a particle is transmitted from left to right of 
the barrier and vice- versa. Similar assignments are done for the scatterer S 2 . 

Thus for the amplitude of transmission we have- t — t l t 2 e lk l e~ al + t 1 r 2 r' 1 t 2 e 3lk l e~ 3al + ... which can be summed 

as 
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FIG. 4: Summing the different paths, Si and S2 denote the two scatterer's. I is the distance between them, e 
denote the propagation and attenuation factors in the locality of interest. 
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and this is the transmission amplitude in presence of wave attenuation. Again for the case of reflection amplitude we 
have r = r 1 + t 1 r 2 i' 1 e 



2ik'l -2al 



+ t 1 r 2 r' 1 r 2 t' 1 e 4ik ' l e- 4al + t 1 r 2 r' 1 r 2 r' 1 r 2 t' 1 e 
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which leads to 



r i ~~ ar 2 e 2lk l e 2al 
r[r 2 t 



(5) 



In Eq. 5, 



1' 1 



t x t\ is the determinant of the S-Matrix of the first barrier and as we are only dealing with 



unitary S-Matrices therefore the determinant is of unit modulus for all barriers. In these expressions k! is the wave 
vector in the region of interest. The transmission and reflection coefficients can be calculated by taking the square of 
the modulus of the expressions in Eq's. (4) and (5). 

The traversal length for transmission and reflection length in case of reflection are calculated as belowQ- 



I = lira n . . 

2a^a d(2a) 

and the reflection length in case of reflection is defined as- 

5(ln|r| 2 ) 



l R = lim -- 



d(2a) 



(6) 



(7) 



The traversal times for transmission or reflection times in case of reflection can be calculated from the formula- 
t R/t _ -y^hgj-gin as before is the speed of propagation in the region of interest. For the case of the potential 

profile sketched in Figure 5, k' — k. From Eq's. (6) and (7) we can calculate the traversal length for transmission- 
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and for reflection- 
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(8) 



(9) 



Here 5ft represents real part of the quantity in brackets. In the above two equations the traversal and reflection lengths 
have been normalized with respect to the length I of the locality of interest, which is the well region of the potential 
profile of Figure 5. Throughout the discussion the quantities are expressed in their dimensionless form. 

We consider the case of two rectangular barriers as shown in Figure 5 separated by a distance I. The energies, 
potentials and the lengths are all in their dimensionless form. 

The S-Matrices for the double heterostructure are given as- 
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FIG. 5: The double heterostructure. Barrier's are denoted by length dj's, height Vj's and j = 1, 2. The length of the well is I. 
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FIG. 6: T and r T for a symmetric double heterostructure . I = 10.0, Vi = V2 = 1.0 and ai = £12 = 0.2. 
wherein K 3 = y/ 2m{ ^r E) , k = ^^and e j± =*-±*i. Herein, j=l,2. 

In Figure 6 we plot the Transmission coefficient and the normalized times t t (=t ) for a symmetrical double 
barrier V = V\ = V2. These times are those required by a quantum particle to traverse or reflect from the well of 
the double heterostructure. We have normalized these times by the time taken by a particle to traverse the length I, 
i.e., St without the potential profile as in Figure 5. From Eq. 8, it is clear that transmission time is always positive. 
Moreover, we can readily show from our treatment that local times spent by the particle in non-overlapping regions 
are additive. As expected the traversal times are larger near the resonances. In case of non-symmetrical structures one 
obtains negative reflection times (which should not be surprising) in some parameter reg ion but traversal times are 
positive for both symmetrical as well as non-symmetrical structures. It can be argued 23] that this is because in case 
of reflection there is a partial wave corresponding to prompt reflection r\ that never samples the region of interest, 
and also this prompt part leads to self interference delays which cause the sojourn time t r to become negative for 
some values of the parameters. If one removes this prompt part, i.e., r np = r — r 1} and we re-calculate the sojourn 
time t r " p with this prompt part removed we find it to be positive and given by r Rnp = t t + 1, as t t is positive. 
In the case of non-symmetrical structures t t is independent of the fact that a particle is incident from the left or 
the right, but t r depends on the direction of the incident particle. There is also a remarkable assertion found in 
the literature|l7j concerning the measurement of the time of transmission or reflection, which is t d = Tt t + Rt r . 

Herein t t and t r are as given above while the dwell time t d = ^ \tp\ 2 dx. tp is the wavefunction in the locality of 
interest and v is the speed of the particle in the region of interest. We have verified explicity after calculation that 
this equivalence does not hold 0, 0] • Our method can also be readily extended to calculation of traversal time of 
tunneling (for case of a single barrier) . Results are in agreement with recent calculations 23J . 
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IV. CONCLUSION 



In conclusion, we have shown that wave attenuation is much better in modeling dephasing due to absorption than 
it's counterpart the optical potential. Further when we extend this method to calculate the conditional sojourn times 
we find here also wave attenuation is easier to deal with than optical potential. 
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